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Very metal-poor stars of the halo are the relics of the star formation in the early Galaxy. Among 
these stars carbon-enhanced metal-poor (CEMP) stars count for 9-25%. Most CEMP stars are 
also enriched in ^-process elements and are called CEMP-s stars. The large number of CEMP- 
s stars observed in binary systems suggests that the chemical enrichment is due to wind mass 
transfer in the past from an AGB donor star on to a low-mass companion, the star that we are 
observing today. However, binary population synthesis models predict CEMP fractions of only 
« 2%. As an alternative to the canonical Bondi-Hoyle-Lyttleton (BHL) wind accretion model, 
recent hydrodynamical simulations suggest an efficient mode of wind mass transfer, called wind 
Roche-lobe overflow (WRLOF), can reproduce observations of AGB winds in binary systems. 
We use our population synthesis model to test the consequences of WRLOF on a population 
of CEMP stars. Compared to previous predictions based on the BHL model we find a modest 
increase of the fraction of CEMP stars and substantial differences in the distributions of carbon 
and periods in the population of CEMP stars. 
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1. Carbon-Enhanced Metal-Poor stars 

The stars that we observe in the Galactic halo are typically low in mass, poor in metals and 
have hardly evolved since their formation. We call very metal-poor (VMP) those stars charac- 
terised by an iron abundance [Fe/H]^ < —2.0. Among the VMP stars of the halo the wide field 
spectroscopic surveys HK [1] and HES [2] observed a large fraction of carbon-enhanced metal-poor 
stars ([C/Fe] > +1.0; CEMP stars). The observed CEMP to VMP ratio varies from 9% to 25% 
[3-5] and increases rapidly for decreasing iron content and increasing distance from the Galactic 
plane [6]. Most observed CEMP stars show an enhancement of heavy elements produced by slow 
neutron-capture process, e.g. barium. CEMP stars with barium abundance [Ba/Fe] > +0.5 are 
called CEMP-5 stai^s and account for approximately 80% of all CEMP stars [7]. 

Three different solutions have been proposed to explain the formation of CEMP stars: (/) 
CEMP stai^s are formed in binary systems where, in the past, the carbon-rich wind of a thermally- 
pulsing asymptotic giant branch (TPAGB) star was accreted to the surface of the companion which 
is the star that we are observing today; (//) CEMP stars produce their own carbon in their interior 
and bring it to the surface by mixing episodes that happen in some stellar evolution models at ex- 
ti'emely low metallicity [8]; (///) the carbon is primordial and was produced by the first generation 
of stars [9]. The binary scenario is the most likely formation mechanism for CEMP-s stars: a spec- 
troscopic analysis shows that the fraction of CEMP-s stars with detected radial-velocity variations 
is consistent with the hypothesis of all being members of binary systems [10]. 

Several studies attempt to reproduce the observed CEMP fraction of 9 — 25% but the typical 
predictions of population synthesis models are at least a factor four lower than the observations 
[11]. Some authors suggest that an initial mass function (IMF) biased towards intermediate-mass 
stars at low metallicity is required to reproduce the CEMP/VMP fraction measured in the halo 
[12, 13], but large changes to the IMF are inconsistent with the small observed fraction of nitrogen- 
enriched metal -poor stars [14]. The IMF at low metallicity is not the only uncertainty that affects 
the study of the origin of CEMP stars; other physical parameters are poorly constrained, e.g. the 
internal mixing and diffusion processes in stars, the nucleosynthesis during the AGE phase, the 
mass-ti"ansfer process, the binary fraction of low-mass stars in the halo and the initial distribution 
of the orbital parameters. For a description of these uncerainties in the context of CEMP stars we 
refer to the work of Izzard and collaborators and references therein [11]. In this work we investigate 
the wind mass-transfer process and its effects on a population of CEMP stars. 

2. Wind mass transfer in binary stars 

In the late stage of its evolution a star of mass between 0.8 — 8M0 becomes a TPAGB star. 
In this phase a sufficiently massive star produces carbon and .s-elements in its interior, brings them 
to the surface by third dredge up and expels them in the strong wind mass loss characteristic of 
TPAGB stars. If the wind-losing star belongs to a binary system, a fraction of the ejected material 
is accreted to the companion stai^ and pollutes its surface. The TPAGB star then evolves into a 



' given two elements X and Y, their abundance ratio is [X/Y] = log j q (Nx /Ny ) — log j q {Nxq /Nyq ) , where Nx j refers 
to the number density of the elements X and Y and denotes the abundance in the Sun. 
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Figure 1: Left panel. Schematic picture of the WRLOF mechanism: the wind acceleration radius R^ lies 
close to the Roche-lobe radius 7?l,i of the AGB star (in red). The wind is slow inside the wind acceleration 
zone (orange area) and can be efficiently accreted to the low-mass secondary through the inner Lagrangian 
point Li (sizes are not in scale). Right panel. Efficiency of wind mass transfer j3acc as a function of the initial 
period for a binary system with Mi = 1 .5 M© and M2 = 0.6 M© . In this and in the plots that follow we show 
models calculated with BHL and WRLOF prescriptions in blue and red respectively. 



white dwarf and today we are able to observe only the chemically-enriched secondary star. In this 
picture the wind mass transfer therefore has a fundamental role for the formation of CEMP stars. 

The canonical description of wind mass transfer is the Bondi-Hoyle-Lyttleton (BHL) pre- 
scription [15], which predicts low accretion efficiencies. The BHL prescription is based on the 
assumption that the wind velocity is much larger than the orbital velocity of the accreting star. 
However this assumption is not always valid in the case of AGB winds. Outflows from AGB stars 
are observed with velocities in the range 5 — 30 kms^^ while the orbital velocity in binary stars of 
periods around 10^ days is about 10kms^\ i.e. comparable. 

Recent hydrodynamical simulations suggest a mode of mass transfer called "wind Roche-lobe 
overflow" (WRLOF) applicable in the case of the slow AGB winds. WRLOF occurs in systems 
where the radius of the wind acceleration zone, i.e. the region where the wind is accelerated beyond 
the escape velocity, is larger than, or is a significant fraction of, the Roche-lobe radius of the wind- 
losing star. When this condition is fulfilled the wind is gravitationally confined to the Roche lobe 
of the donor star and is funnelled through the inner Lagrangian point L\ towards the secondary, on 
to which it accretes with high efficiency (up to 50%) over a wide range of separations [16]. We 
schematically represent the conditions in which the WRLOF mode of mass transfer occurs in the 
left panel of Fig. 1 . 

AGB winds are thought to be accelerated by a combination of stellar pulsation and radiation 
pressure. The pulsations increase the gas density in the cold outer atmosphere where dust grains 
form. The dust grains are accelerated by the radiation pressure and drag the surrounding gas along 
due to coUisional momentum transfer [17]. Therefore, we assume that the wind acceleration radius 
coincides with the dust formation radius R^, which depends on the stellar radius, on the effective 
temperature of the star and on the condensation temperature of the dust (the higher this temperature 
is, the closer to the stellar surface the dust forms). During most of the AGB phase the dust formation 
radius is linearly proportional to the stellar radius and for carbon-rich dust R^/R^ ~ 3. 
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In this work we implement in our binary population synthesis code a model of the WRLOF 
process based on the results of the above-mentioned hydrodynamical simulations [18]. We present 
the results of the comparison between the evolution of a population of binary systems in the case 
of wind mass transfer modelled by the BHL prescription and by WRLOF. 

3. Results 

Our population synthesis simulations consist of A^^ binary stars in InMi — lnM2 — \na pa- 
rameter space, where Mi 2 are the initial masses of the primary and the secondary star, a is the 
initial separation and we choose A'^ = 128. The initial masses M\ and M2 vary respectively in the 
range [0.8, 8.0] M0 and [0.1, O.QJMq. The initial separation is chosen in the range [3, 10^] R© and 
we consider a binary fraction of unity in this range. The metallicity is Z = 10^^ which corre- 
sponds to [Fe/H] Ri —2.3. Stars are selected from our model population according to their age, 
surface gravity and surface abundances: VMP stars are older than 10 Gyr and with surface gravity 
logio(g/cms^^) < 4.0, i.e. we select essentially subgiants or giants; CEMP stars are VMP stars 
with [C/Fe] > +1.0. For a complete discussion of our model, its characteristics, parameters and 
uncertainties we refer to the work of Izzard and collaborators and our forthcoming paper [11, 18]. 

In the right panel of Fig. 1 we show an example of the effect of our model of WRLOF 
on the accretion efficiency of wind mass transfer in a binary system with M\ = 1.5 M0 , M2 = 
O.6M0 and initial period in the range [10^ , 10^] days. We compare the accretion efficiency j8acc 
calculated with the BHL prescription (blue line) and with our WRLOF model (red line). The 
maximum accretion efficiency in the case of WRLOF is more than 40% and a high efficiency is 
possible over a wide range of periods. A 1.5M0 stai^ loses approximately O.8M0 during the AGE 
phase, therefore its 0.6 M0 companion can accrete enough mass to become carbon-enhanced. The 
maximum accretion efficiency in the BHL model is about 10% for periods around 3000 days, which 
means the companion accretes 0.08 M0 or less. A 0.68 M0 star is not massive enough to evolve off 
the main sequence and is not selected as a CEMP star by our model. 

The effect of WRLOF on a population of binary stars is to widen towards longer periods and 
lower secondary masses the range of initial separations and masses of systems that become CEMP 
stars. As a result we calculate a CEMP A'^MP ratio of 2.63% —4.06%, an increase by a factor 1.2 — 
1.8 compared to prior studies [11], where the range of variations is due to the assumptions made on 
the amount of angular momentum carried away by the wind lost and on the dependence of the wind 
accretion efficiency on the mass ratio of the two stars. The CEMP fraction is still small compared 
to the observed fraction of CEMP stars. In Fig. 2 (left panel) we show the distribution of the 
carbon abundance calculated with the BHL (blue) and the WRLOF (red) prescriptions compared 
with the observations (histograms with Poisson en^ors). Our observed sample is based mainly on 
the SAGA observational database and contains all data on CEMP stars currently available in the 
hterature [3, 10, 19]. The WRLOF model shows qualitatively the same trend as the data, but 
peaks at a carbon abundance higher than the observations. This may be an indication that in our 
model the accretion efficiency is too high. The transferred material is mixed in the envelope of 
the secondary star, composed mainly of hydrogen. A star of 0.4 M0 must accrete at least other 
0.4 M0 from the donor to be selected as a CEMP star and this material is only weakly diluted in 
its envelope. Hence secondary stars with low initial mass become strongly carbon-emiched stars. 
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Figure 2: Distribution of [C/Fe] {left panel) and periods {right panel) in the CEMP population. Histograms 
show the observed distribution in our data sample, with Poisson errors. Plus signs are periods of observed 
CEMP stars. The bins are equally spaced and the width of each bin is 0.2 in [C/Fe] and logiQ(P/days) 
respectively. The y-axis indicates the expected number of CEMP stars in each bin. The plot is normalised 
such that the area under the graph is the same and BHL model peaks at 1 . 

If we reduce the accretion efficiency of tlie mass transfer very low-mass stars do not increase tiieir 
mass enougli to be selected as CEMP stars and tliis results in a shift of the distribution towards 
lower values of [C/Fe]. In Fig. 2 we compare the final period distributions predicted by our models 
with the observations. Even though many CEMP stars are known to be in binaries, only a few have 
known orbital periods because long orbital periods are difficult to measure. Our model based on 
the WRLOF prescription predicts most CEMP stars at orbital periods around 2000 days, more than 
one order of magnitude less than the model based on the BHL prescription. Both our models have 
difficulty reproducing the shortest-period CEMP stars. 

4. Summary 

We include for the first time the WRLOF process of mass transfer in a binary population syn- 
thesis model and we study the effects on a population of CEMP stars. Our WRLOF model predicts 
an increase of the CEMPA'^MP ratio up to 4.06%, a factor 1.8 higher than the predictions with the 
BHL prescription. However, this ratio is still low compared to the observed CEMP fraction. The 
WRLOF period distribution peaks at shorter periods compared to the canonical BHL model but 
the two shortest-period CEMP stars are not reproduced. The carbon distribution predicted by our 
WRLOF model is qualitatively similar to the observations but it peaks at a higher value of [C/Fe]. 
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